Search for photon line-like signatures from Dark Matter annihilations with H.E.S.S. 
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Gamma-ray line signatures can be expected in the very-high-energy (VHE; E 1 > 100 GeV) do- 
main due to self-annihilation or decay of dark matter (DM) particles in space. Such a signal would 
be readily distinguishable from astrophysical 7-ray sources that in most cases produce continu- 
ous spectra which span over several orders of magnitude in energy. Using data collected with the 
H.E.S.S. 7-ray instrument, upper limits on line-like emission are obtained in the energy range be- 
tween ~ 500 GeV and ~ 25 TeV for the central part of the Milky Way halo and for extragalactic 
observations, complementing recent limits obtained with the Fermi-LAT instrument at lower ener- 
gies. No statistically significant signal could be found. For monochromatic 7-ray line emission, flux 
limits of (2 x 10~ 7 - 2 x 10" 5 ) m" 2 s~ 1 sr" 1 and (1 x 10" 8 - 2 x 10~ 6 ) m" 2 s" 1 sr" 1 are obtained for 
the central part of the Milky Way halo and extragalactic observations, respectively. For a DM par- 
ticle mass of 1 TeV, limits on the velocity-averaged DM annihilation cross section {av) xx ^ 11 reach 
cm 3 s _1 , based on the Einasto parametrization of the Galactic DM halo density profile. 



INTRODUCTION 



In the last few years, imaging atmospheric Cherenkov 
telescopes (IACTs) have been used to search for dark 
matter (DM) signals in very-high-energy (VHE; E 1 > 
100 GeV) 7 ravs jlWIojl. Objects with large predicted DM 
density, like the Galactic centre (GC), the central Galac- 
tic halo region (CGH), dwarf galaxies or centres of nearby 
galaxies were studied. All such searches concentrated on 
the detection of 7 rays produced in decays of secondary 
particles - mostly neutral mesons - in the process of DM 
self-annihilation or decay (see, e.g., [Uj, • The broad 
energy distribution of such 7 rays is continuous and there- 
fore more difficult to distinguish from 7-ray emission from 
astrophysical (particle accelerating) sources, as opposed 
to spectral features, which would pose a much more strik- 
ing evidence for a DM-induced 7-ray signal. The most 



prominent spectral feature is a 7-ray lineQ, which, for 
DM self-annihilation into 77/7.Z (and m x 3> TOz), is ex- 
pected at an energy at or close to the DM particle mass, 
E 1 k. m x . For a decay x 1% of a DM particle x with 
m x 3> mx, Ey w m x /2. Such annihilations or decays 
are, however, loop-suppressed, since electrically neutral 
DM particles do not couple to photons directly. Nonethe- 
less, recent theoretical developments show the possibil- 
ity of a rather pronounced spectral feature for some im- 
plementations of particle physics beyond the Standard 
Model (see, e.g., [FJ|). In these models, spectral signa- 
tures may arise due to the radiation of a hard photon 
from real or virtual charged particles created in the an- 
nihilation process and be a dominant component in the 



Note, however, that VHE 7-ray line features may also arise due 
to unshocked e+/e - -winds created by pulsars [13|| . 



3 



overall 7-ray annihilation spectrum. Here a search for 7- 
ray line-like signatures conducted with the H.E.S.S. ex- 
periment in the energy range E y ~ 500 GeV — 25 TeV 
is reported, complementing a recent search at energies 
between 7 GeV and 200 GeV with the Fermi-LAT instru- 
ment 15J and studies discussing an indication for a line 
feature at an energy of about 130 GeV [IMIj]. 

The search for a DM-induced spectral signature in the 
H.E.S.S. data is performed separately for two sky regions 
of interest. The first is the CGH, a promising region due 
to its proximity and predicted large DM concentration. 
Following jj§j, the search region is defined as a circle of 
1° radius centred on the GC, where the Galactic plane is 
excluded, by requiring |6| > 0.3°. The second region is 
the extragalactic sky covered by H.E.S.S. observations, 
with regions containing known VHE 7-ray sources being 
excluded from the analysis. For both data sets, the un- 
certainty on the strength of a putative DM annihilation 
signal is much reduced in comparison to the observations 
of centres of galaxies: for the CGH, the very centre is not 
considered, thus avoiding a region where the DM profile 
is only poorly constrained [§| . For the extragalactic data 
set, differences in DM density between individual sub- 
structures are averaged out by observing many different 
fields of view {TgJ. One should note, however, that a 
potentially large (but highly uncertain) 7-ray flux from 
Galactic DM annihilations may contribute to the extra- 
galactic analysis [2(3]. 
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FIG. 1. Reconstructed flux spectrum of the CGH region, us- 
ing 25 equidistant bins per unit of log 10 (_E 7 ). Flux points have 
been multiplied by E^' 1 . The data consist mostly of hadronic 
cosmic ray background events, reconstructed using a 7-ray hy- 
pothesis. The spectrum is well described by the parametriza- 
tion introduced in Eq. [TJ depicted by the black solid line. The 
corresponding x 2 -test probability is p = 0.34. The two contri- 
butions P(x) and G(x) are shown by the dashed-dotted and 
the dashed curve, respectively. Note that the shape of the 
Gaussian function G(x) is much broader than the expected 
monochromatic line feature from DM annihilations. As an 
example, the red curve shows the expected signal of a line at 
E-, = 2 TeV that would be detected with a statistical signifi- 
cance of 5 standard deviations above the background. 



METHODOLOGY AND RESULTS 

The CGH data set is composed of 112 h (live time) of 
GC observations recorded with the H.E.S.S. VHE 7-ray 
instrument (see J2l| and references therein) during the 
years 2004-200qj. The mean distance between the tele- 
scope pointing positions and the GC is 0.7°, with a max- 
imum of 1.5° [8|. The extragalactic data set comprises 
1153 h of H.E.S.S. observations taken during 2004-2007, 
targeted at various extragalactic objects. Regions in the 
field-of-view (FoV) containing known VHE 7-ray sources 
are excluded by masking out a circular region (of radius 
0.2° for point sources) around the source position. 

Observations with zenith angles larger than 30° are 
excluded from the analysis to lower the energy thresh- 
old, resulting in a mean zenith angle of 14° (19°) for the 
CGH (extragalactic) observations. Only 7-ray-like events 
are accepted for which the distance between the recon- 
structed 7-ray direction and the observation direction of 
the H.E.S.S. array is smaller than 2°, avoiding showers 



2 Data from later periods were excluded, since the gradual degra- 
dation in time of the optical efficiency of the instrument would 
result in an increased energy threshold. 



being reconstructed too close to the edges of the ~ 5° di- 



ameter FoV of the H.E.S.S. cameras [2lJ. Furthermore, 
events are considered only if they pass H.E.S.S. standard 
7-ray selection criteria defined in [2l| and triggered all 
four telescopes. Only 15 % of the total event sample is 
kept by the latter selection. However, compared to the 
H.E.S.S. standard analysis, such selection leads to a bet- 
ter signal to background ratio and an improved energy 
resolution of Gaussian width cte (17% at 500 GeV and 
11% at 10 TeV), and therefore increases the sensitivity 
of the analysis to spectral features by up to 50%. The 
energy threshold is 310 GeV (500 GeV) for the CGH (the 
extragalactic) data set. 

Differential flux spectra are calculated from the re- 
constructed event energies separately for the CGH and 
extragalactic data sets using zenith angle-, energy- and 
offset-dependent effective collection areas from 7-ray sim- 
ulations. Since sky regions containing known VHE 7-ray 
sources were excluded from the analysis, the spectra con- 
sist mostly of 7-ray-like cosmic-ray background events 
(and a fraction of ~ 10% of electrons). These spectra are 
well described by the empirical parametrization 
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where E 7 is the reconstructed energy of the event under 
7-ray hypothesis and P(x) — exp (a\x + a,2X 2 + 03a; 3 ). 
G{x) is a Gaussian function with mean /i x and rms ct x , 
and x = log 10 (-E 7 /l TeV). The free parameters a ...3, j3, 
/i x , and ct x are optimized simultaneously by a maximum 
likelihood approach based on the binned event count 
spectrum. Since the number of reconstructed counts rii in 
energy bin i of the count spectrum is Poisson-distributed, 
the log-likelihood function takes the form 

N 



In C = rii In Aj — Aj 



where A^ is the number of counts in bin i that is expected 
according to the flux spectrum parametrization given in 
Eq. [TJ and N is the total number of bins of the count 
spectrum. As an example, Fig. [T] shows the differential 
flux spectrum and the best-fit background parametriza- 
tion obtained for the CGH data set. 

On top of the smooth cosmic ray flux spectrum, a 
monochromatic 7-ray lin^f] may be identified as a Gaus- 
sian peak of width cte centred at the line energy E 7 . To 
search for such lines, a Gaussian term with fixed energy 
E-y and fixed corresponding width cte was added to the 
spectrum parametrization given in Eq. [T] The spectrum 
was refit, and from the normalization of the Gaussian the 
flux of the putative line was reconstructed. By repeat- 
ing this procedure, using ten logarithmically equidistant 
energies E~ l per decade of energy, the flux spectrum was 
scanned for monochromatic 7-ray signatures. Line scans 
were performed in the energy range 0.5TeV-20TeV and 
0.8 TeV-25 TeV for the CGH and the extragalactic data 
sets, respectively. 

No 7-ray line flux was found to exceed the a-priori 
chosen detection threshold of A In C = 12.5, correspond- 
ing to a significance of 5 standard deviations above the 
background level for Gaussian parameters. Thus flux up- 
per limits were calculated by constraining the flux nor- 
malization of the Gaussian to be non-negative in the fit 
and using the MINOS package from the Minuit[22| fit- 
ting tool to calculate asymmetric errors with error level 
Aln£ = 1.35, corresponding to a 95% CL one-sided limit 
on the flux of the line 15|, l23[ . These limits are shown 
in Fig. [2j To test whether the limits are compatible with 
random fluctuations of the background, a large number 
of statistically randomized fake background spectra was 
simulated using the best-fit background parametrization 
as an input, and limits were obtained for each of these 
spectra. The resulting mean limits, together with the 
68% CL region calculated from the limit distribution at 
each test energy, are shown in Fig.[5]for comparison. Also 



shown are mean reconstructed fluxes from simulated lines 
that are detected with a significance of 5 standard devi- 
ations using the above prescription. 

Additionally, flux upper limits were determined for 
broader spectral features like those arising due to inter- 
nal bremsstrahlung (IB). As an example, calculations by 



14] in the framework of supersymmetric models predict 



the contribution of IB photons to the 7-ray spectrum 
to dominate over secondary 7-ray production for photon 
energies close to the DM (neutralino) mass m x . Flux 
up per limits for the benchmark models BM2 and BM4 
of [l4j were calculated following the technique described 
above. Firstly, the signal shapes predicted by the models 
were convolved with the energy response of the instru- 
ment. Together with the background parametrization, 
the resulting templates were then fitted (with the nor- 
malization of the template and the background parame- 
ters being free variables in the fit) to the flux spectrum. 
Note that only the IB part of the full annihilation spec- 
tra of these models is considered since the contribution 
from production of secondary photons steeply decreases 



towards m x (see |14[), and is therefore hard to discrim- 



3 In this context, the term 'monochromatic line' refers to spec- 
tral features with energy width much smaller than the energy 
resolution eg of the H.E.S.S. instrument. 



inate against the cosmic-ray background. In any case, 
since these models were calculated for a very specific set 
of MSSM parameters (and hence neutralino mass), they 
can only serve as a template to demonstrate the sensi- 
tivity of H.E.S.S. for features of similar shape (and are 
therefore referred to as BM2-like and BM4-like limits). 
Fig. [3] shows that - because of the intrinsic widths of the 
expected features - these limits are typically weaker by 
a factor two (BM2-likc) to ten (BM4-like) compared to 
the monochromatic line limits. Note that all flux lim- 
its do also constrain putative features in the spectrum of 
cosmic ray electrons and positrons, since the H.E.S.S. ex- 
periment exhibits a similar sensitivity for detecting these 
particles as for 7 rays. 

Possible systematic uncertainties due to the unknown 
shape of the background spectrum have been extensively 
studied, e. g. by changing the background parametriza- 
tion described in Eq. [1] to one based on Legendre poly- 
nomials. The background parametrization does not show 
any significant correlation with shape parameters of spec- 
tral signatures, in particular with regard to the G(x) 
term. The stability of the 7-ray flux reconstruction was 
investigated by adding artificial peaks to the background 
spectrum and reconstructing them with the fitting pro- 
cedure described above. The systematic uncertainty on 
the reconstructed peak flux was of the order of a few 
percent, and the fit of the background was found to be 
very stable and independent of the location and normal- 
ization of the artificial peak. On the other hand, despite 
detailed Monte-Carlo simulations of the instrument, the 
true energy resolution cte of the instrument might be un- 
derestimated. When cte is artificially enlarged by e. g. 
20% - i.e. cte = 20% (13%) at E 1 = 500 GeV (10 TeV) 
-, upper limits get shifted to larger values by about 15- 
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FIG. 2. Upper limits on 7-ray flux from monochromatic line 
signatures, derived from the CGH region (red arrows with 
full data points) and from extragalactic observations (black 
arrows with open data points). For both data sets, the solid 
black lines show the mean expected limits derived from a large 
number of statistically randomized simulations of fake back- 
ground spectra, and the gray bands denote the corresponding 
68% CL regions for these limits. Black crosses denote the flux 
levels needed for a statistically significant line detection in the 
CGH dataset. 
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FIG. 3. Flux upper limits on spectral features arising from 
the emission of a hard photon in the DM annihilation pro- 
cess. Limits are exemplary shown for features of comparable 
shape to those arising in the models BM2 and BM4 given in 



[14j. The monochromatic line limits, assuming m x 
shown for comparison. 



E~, are 



20%, depending on the energy and the statistics in the 
individual spectrum bins. The maximum shift is ob- 
served in the extragalactic limit curve and amounts to 
40%. In total, the systematic error on the flux upper 
limits is estimated to be about 50 %. All flux upper 
limits were cross-checked using an alternative analysis 
framework 2J] , with an independent calibration of cam- 
era pixel amplitudes, and a different event reconstruction 



FIG. 4. Limits on the velocity-weighted cross section for DM 
annihilation into two photons calculated from the CGH flux 
limits (red arrows with full data points). The Einasto density 
profile with parameters described in [2(j] was used. Limits ob- 
tained by Fermi-LAT, assuming the Einasto profile as well, are 
shown for comparison (black arrows with open data points) 

M. 



and event selection method, leading to results well con- 
sistent within the quoted systematic error. 

For the Einasto parametrization of the DM density 
distribution in the Galactic halo [2oj |. limits on the 
velocity-weighted DM annihilation cross section into 7 



rays, (av) > 



are calculated from the CGH flux limits 



using the astrophysical factors given in [8( . The result is 
shown in Fig. [4] and compared to recent results obtained 
at GeV energies with the Fermi-LAT instrument. 



SUMMARY AND CONCLUSIONS 

For the first time, a search for spectral 7-ray signatures 
at very-high energies was performed based on H.E.S.S. 
observations of the central Milky Way halo region and ex- 
tragalactic sky. Both regions of interest exhibit a reduced 
dependency of the putative DM annihilation flux on the 
actual DM density profile. Upper limits on monochro- 
matic 7-ray line signatures were determined for the first 
time for energies between ~ 500 GeV and ~ 25 TeV, cov- 
ering an important region of the mass range of particle 
DM. Additionally, limits were obtained on spectral sig- 
natures arising from internal bremsstrahlung processes, 
as predicted by the models BM2 and BM4 of It 
should be stressed that the latter results are valid for 
all spectral signatures of comparable shape. Besides, all 
limits also apply for potential signatures in the spectrum 
of cosmic-ray electrons and positrons. 

Flux limits on monochromatic line emission from the 
central Milky Way halo were used to calculate upper lim- 
Limits are obtained in a neutralino 



its on (av) 
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mass range that is complementary to the line searches 
performed by Fermi-LAT [TBI, reaching ~ 10~ 27 cm 3 s -1 
at a DM mass of 1 TeV, comparable to those obtained by 
Fcrmi-LAT at energies of ~ 100 GeV. 
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